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Unequal Crossover Dynamics 
in Discrete and Continuous Time 

Abstract. We analyze a class of models for unequal crossover (UC) of sequences con- 
taining sections with repeated units that may differ in length. In these, the probability of 
an 'imperfect' alignment, in which the shorter sequence has d units without a partner in 
the longer one, scales like q d as compared to 'perfect' alignments where all these copies are 
paired. The class is parameterized by this penalty factor q. An effectively infinite population 
size and thus deterministic dynamics is assumed. For the extreme cases q — and 5=1, 
and any initial distribution whose moments satisfy certain conditions, we prove the conver- 
gence to one of the known fixed points, uniquely determined by the mean copy number, in 
both discrete and continuous time. For the intermediate parameter values, the existence of 
fixed points is shown. 



1. Introduction 

Recombination is a by-product of (sexual) reproduction that leads to the mixing 
of parental genes by exchanging genes (or sequence parts) between homologous 
chromosomes (or DNA strands). This is achieved through an alignment of the 
corresponding sequences, along with crossover events which lead to reciprocal ex- 
change of the induced segments. In this process, imperfect alignment may result 
in sequences that differ in length form the parental ones; this is known as unequal 
crossover (UC). Imperfect alignment is facilitated by the presence of repeated el- 
ements (as is observed in some rDNA sequences, compare Q), and is believed 
to be an important driving mechanism for their evolution. The repeated elements 
may follow an evolutionary course independent of each other and thus give rise to 
evolutionary innovation. For a detailed discussion of these topics, see (22 23 1 and 
references therein. 

This article is concerned with a class of models for UC, originally investigated 
by Shpak and Atteson 1 22 1 for discrete time, which is built on preceding work 
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Fig. 1. Snapshot of an unequal crossover event as described in the text. Rectangles denote 
the relevant blocks, while the dashed lines indicate possible extensions with other elements 
that are disregarded here. 



by Ohta [ 15 1 and Walsh [25 1 (see [22 1 for further references). Starting from their 
partly heuristic results, we prove various existence and uniqueness theorems and 
analyze the convergence properties, both in discrete and in continuous time. This 
will require a rather careful mathematical development because the dynamical sys- 
tems are infinite dimensional. 

In this model class, one considers individuals whose genetic sequences contain 
a section with repeated units. These may vary in number, i E Nq — {0, 1,2,3,...}, 
where i = is explicitly allowed, corresponding to no unit being present (yet). The 
composition of these sections (with respect to mutations that might have occurred) 
and the rest of the sequence are ignored here. 

In the course of time, recombination events happen, each of which basically 
consists of three steps. First, independent pairs are formed at random (in equidis- 
tant time steps, or at a fixed rate). Then, their respective sections are randomly 
aligned, possibly with imperfections in form of 'overhangs', according to some 
probability distribution for the various possibilities. Finally, both sequences are 
cut at an arbitrary common position between two adjacent building blocks, with 
uniform distribution for the cut positions, and their right (or left) fragments are in- 
terchanged. This so-called unequal crossover is schematically depicted in Figure^ 
Obviously, the total number of relevant units is conserved in each event. 

We assume the population size to be (effectively) infinite. (Concerning finite 
populations, see the remarks in Section0) Then, almost surely in the probabilistic 
sense, compare |4, Sec. 11.2], the population is described by the deterministic 
time evolution of a probability measure p S M.^ (No), which we identify with 
an element p = (j5/ c ) J tgN m me appropriate subset of £ 1 (No). Since we will not 
consider any genotype space other than No in this article, reference to it will be 
omitted in what follows. These spaces are complete in the metric derived from the 
usual i 1 norm, which is the same as the total variation norm here. The metric is 
denoted by 




(1) 



With this notation, the above process is described by the recombinator 
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Here, TJy ^ denotes the probability that, given a pair (k, £), this pair turns into 
(i, j). Consequently, for normalization, we require 

£)ly,M = l forallfc^eNo. (3) 

The factor in (0 describes the probability that a pair (fc,£) is formed, i.e., 
we assume that two individuals are chosen independently from the population. We 
assume further that ^ = 2},-^ = Ty.^-, i.e., that T^./^ is symmetric with 
respect to both index pairs, which is reasonable. Then, the summation over j in Q 
represents the breaking-up of the pairs after the recombination event. These two 
ingredients of the dynamics constitute what is known as (instant) mixing and are 
responsible for the quadratic nature of the iteration process. 

As mentioned above, we will only consider processes that conserve the total 
copy number in each event, i.e., T^ M > for i + j = k + £ only. Together 
with the normalization Q and the symmetry condition from above, this yields the 
(otherwise weaker) condition 

v - "* • rp i+j k + £ 

which implies conservation of the mean copy number in the population, 

Condition Q and the presence of the pref actor l/||p||i in (13 make 1Z norm 
non-increasing, i.e., HT^a^Hj ^ ||aj||i, and positive homogeneous of degree 1, 
i.e., IZ(ax) = \a\lZ(x), for all x € £ l and a 6 ML Furthermore, 1Z is a positive 
operator with ^(xjHj = ||x||j for all positive elements x G t 1 . Thus, it is 
guaranteed that 1Z maps M.f> the space of positive measures of total mass r, into 
itself. This space is complete in the topology induced by the norm i.e., by 
the metric d from Q. (For r = 1, of course, the prefactor is redundant but ensures 
numerical stability of an iteration with 7Z.) 

Given an initial configuration p Q = p(t = 0), the dynamics may be taken in 
discrete time steps, with subsequent generations, 

p{t + l)=K{p(t)), teN . (5) 

Our treatment of this case will be set up in a way that also allows for a generaliza- 
tion of the results to the analogous process in continuous time, where generations 
are overlapping, 

f t p(t) = g(K-l)(p(t)), i£% . (6) 

Obviously, the (positive) parameter q in l|6} only leads to a rescaling of the time t. 
We therefore choose q = 1 without loss of generality. Furthermore, the formula- 
tion of discrete versus continuous time dynamics in Q and (|6j is chosen so that 
the fixed points of Q are identical to the equilibria of regardless of q. This 
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can easily be verified by a direct calculation. In what follows, we will thus use the 
term fixed point for both discrete and continuous dynamics. 

In the UC model, one distinguishes 'perfect' alignments, in which each unit in 
the shorter sequence has a partner in the longer sequence, and 'imperfect' align- 
ments, with 'overhangs' of the shorter sequence relative to the longer one. To come 
to a reasonable probability distribution for the various possiblities, the first are 
taken to be equally probable among each other, while the latter are penalized by 
a factor q d relative to the first, where q £ [0, 1] is a model parameter and d is 
the length of the overhang (at most the entire length of the shorter sequence; in 
the example of Figure ^ interpreted as a snapshot right after the crossover event 
took place, we have d = 1). In the extreme case q = 0, only perfect alignments 
may occur, whereas for q = 1 overhangs are not penalized at all and one obtains 
the uniform distribution on the possibilities. For obvious reasons, the first case is 
dubbed internal UC, the second random UC \ 22\. 

It is now straightforward, though a bit tedious, to derive the transition proba- 
bilities T^j\ e . To this end, one has to trace what happens in steps two and three of 
the recombination event only, while the random formation of pairs does not enter 
here. This has been done in |22| and need not be repeated. However, in view of 
our above remarks, it is desirable to rewrite the findings in a way that reflects the 
natural symmetry properties of the T^\. e . In compact notation, this leads to the 
transition probabilities 

T^u = S i+jtk+i (1 + min{fc, £,i,j}) g«M*Ai-iAJ) , (7) 

where k V £ := max{fc, £}, k A I := min{fc, £}, and 0° = 1. The Cjf) are chosen 

such that holds, i.e., J2i j>o ^ifht = •"■> an( ^ me nence symmetric in k and I. 
Explicitly, they read (see also |22 Sec. 2.1]) 

= (1-g) 2 

ke (k A £ + l)(\k -£\ + 1)(1 - qf + 2q{k A £ - (k A £ + l)q + qk^+i) ' 

Note further that the total number of units is indeed conserved in each event and 
that the process is symmetric within both pairs. Hence ® is satisfied. 

Let us briefly come back to the question of 'discrete' versus 'continuous' time, 
which are considered simultaneously for good reasons. Common to both is the 
nonlinearity that stems from the probability that a certain (random) pair is formed 
in the first place. Then, for the discrete time dynamics Q, the T^\ t have the direct 
meaning of the transition probability that, given a pair (k, £), this turns into a pair 
(i, j). In contrast, for the continuous time dynamics (|6), the number T^ M is to be 
understood as the probability to obtain a pair conditioned on a recombination 
event with a pair of type (k,£), of which each recombines at the same rate. In 
probabilistic terminology, the 1Z of (|5} is the discrete time skeleton of the process 
in continuous time, also called the embedded discrete time process. 

The aim of this article is to find answers to the following questions: 
1. Are there fixed points of the dynamics? 
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2. Given the mean copy number m, is there a unique fixed point? 

3. If so, under which conditions and in which sense does an initial distribution 
converge to this fixed point under time evolution? 

Of course, the trivial fixed point with p = 1 and p k = for k > always exists, 
which we generally exclude from our considerations. But even then, the answer to 
the first question is positive for general operators of the form Q that satisfy Q 
and some rather natural further condition. This is discussed in Section|2] For the 
extreme cases q — (internal UC) and q — 1 (random UC), fixed points are known 
explicitly for every in and it has been conjectured 1 22 1 that, under mild conditions, 
also questions 2 and 3 can be answered positively for all values of q 6 [0, 1]. 
Indeed, for both extreme cases, norm convergence of the population distribution 
to the fixed points can be shown, which is done in Sections[3]and[5l respectively. 
Since the dynamical systems involved are infinite dimensional, a careful analysis 
of compactness properties is needed for rigorous answers. The proofs for q = 1 
are based on alternative representations of probability measures via generating 
functions, presented in Section 0] For the intermediate parameter regime, we can 
only show that there exists a fixed point for every m, but neither its uniqueness nor 
convergence to it, see Section|6] Some remarks in Section[7]conclude this article. 

2. Existence of fixed points 

Let us begin by stating the following general fact. 

Proposition 1. If the recombinator 1Z of satisfies (|3j, then the global Lipschitz 
condition 



is satisfied, with constant C = 3 on £ , respectively C = 2 if x, y € A4 r . 
Proof. Let x, y £ I 1 be non-zero (otherwise the statement is trivial). Then, 



\\'R(x)-'R{y)\\ 1 ^C\\x-y\ 




I|a5-y|li + 



1 



llz/lli^- Nli2/|| r 




The last term becomes 
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7r-jj-||l|y|li^ - IMIii/lli 
\\ x \\i 

= -iT^n — II (II 2/11 1 - \\ x \\i) x + \\ x \\i( x - y)||i ^ 2 ll*-i/lli, 
ii x ii i 

from which ||72.(aj) — lZ(y)\\ i ^ 3||cc — y\\i follows for x, y G I 1 . If a;, y € Ai r , 
the above calculation simplifies to ||7£(a;) — 7^(^)1^ ^2||a; — y|| x . □ 

In continuous time, this is a sufficient condition for the existence and uniqueness 
of a solution of the initial value problem (|5Jl, cf. 1 1 Thms. 7.6 and 10.3]. Another 
useful notion in this respect is the following. 

Definition 1. fl] Sec. 18] Let Y be an open subset of a Banach space E and let 
f : Y — > E satisfy a (local) Lipschitz condition. A continuous function L from 
X C Y to Wis called a Lyapunov function for the initial value problem 

f t x(t) = f(x(t)), x(0) = x eX, 

if the orbital derivative L(xq) ■= liminf t ^ + \{L{x(t)) — L(cco)) satisfies 

L(x ) < (8) 

for all initial conditions Xq € X. 

If further L(xf) — for a single fixed point xy only, then L is called a strict 
Lyapunov function. If a Lyapunov function exists, we have 

Theorem 1. [l] Thm. 17.2 and Cor. 18.4] With the notation of Definition^ assume 
that there is a Lyapunov function L, that the set X is closed, and that, for an initial 
condition Xq G X, the set {x(t) : t S W^a,x(£) exists} is relatively compact in 
X. Then, x{t) exists for allt^O and 

lim dist(x(t),X L ) = 0, 

t — >oo 

where dist(a;,Xi) = mf ye x L \\ x — y\\ and Xl denotes the largest invariant 
subset of {x G X : L{x) = 0} (in forward and backward time). □ 

Obviously, if L is a strict Lyapunov function, we have Xl = { x f} and this theo- 
rem implies d(x(t), xp) — > as t — > oo. 

Returning to the original question of the existence of fixed points, we now 
recall the following facts, compare |3 21 1 for details. 

Proposition 2. 1271 Cor. to Thm. V.1.5] Assume the sequence (p^ n n in A4f to 
converge in the weak-* topology (i.e., pointwise, or vaguely) to some p € M.\, 
i.e., 

lim pi — p k for all k G No, with p k and ^2h>nPh = !• 

77te«, !/ a/so converges weafcfy (in fne probabilistic sense) and in total variation, 
t&,lim n _ o||pW-p|| 1 =a □ 
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Proposition 3. Assume that the recombinator TZ from (|2j satisfies Q and has a 
convex, weak-* closed invariant set M C A4^, i.e., TZ(M) C M, that is tight, i.e., 
for every e > there is an m G No such that X)fc>m Pfc ^ e / or every p S M. 
77ien, 7?. /las a fixed point in M. 

Proof. Prohorov's theorem |21 Thm. III.2.1] states that tightness and relative 
compactness in the weak-* topology are equivalent (see also [3 Chs. 1.1 and 1.5]). 
In our case, M is tight and weak-* closed, therefore, due to Proposition |2] norm 
compact. Furthermore, M is convex and TZ is (norm) continuous by Proposition^ 
Thus, the claim follows from the Leray-Schauder-Tychonov fixed point theorem 
itTslThm. V.191. □ 

With respect to the UC model, we will see that such compact invariant subsets 
indeed exist. 



3. Internal unequal crossover 

After these preliminaries, let us begin with the case of internal UC with perfect 
alignment only, i.e., q = in Q. This case is the simplest because, in each recom- 
bination event, no sequences exceeding the longer of the participating sequences 
can be formed. Here, on A4^ , the recombinator (|2jl simplifies to 

From now on, we write lZ q rather than 1Z whenever we look at a recombinator 
with (fixed) parameter q. It is instructive to generalize the notion of reversibility 
(or detailed balance, compare 1221 (4.1)]). 

Definition 2. We call a probability measure p S A^J" reversible for a recombina- 
tor TZ of the form (|2j if, for all i, j, k, I 0, 

r ^ij,klPkPl =r ^kl,ijViPj- (10) 

The relevance of this concept is evident from the following property. 
Lemma 1. If p G A4 f is reversible for TZ, it is also a fixed point of TZ. 
Proof. Assume p to be reversible. Then, by Q, 

K{p)i= T ij,kiPkPi= T ki, l3 PiP ] =P l ^P ] =P l - 

So, in our search for fixed points, we start by looking for solutions of dl0> . Since, 
for q — 0, forward and backward transition probabilities are simultaneously non- 
zero only if {i,j} = {k,£} C {n,n + 1} for some n, the components p k may 
only be positive on this small set as well. By the following proposition, this indeed 
characterizes all fixed points for q = 0. 
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Proposition 4. A probability measure p G A4J is a fixed point of TZq if and 
only if its mean copy number m = X]fe>o k Pk is finite, P\ m \ = \jn\ + 1 — m, 
P\m] = m + 1 — I" 771 !' an d Pk = for all other k. This includes the case that m 
is integer and p lml = p H = p m = 1. 

Proof. The 'if part was stated in 1221 Sec. 4.1] and follows easily by insertion into 
or JlOb . For the 'only if part, let i denote the smallest integer such that p i > 0. 
Then, 

ft(p)i = Pi + 2 Pi T+i ^ Pl \ Pl + Pl+1 + ^ ITT Pi+e ^ Pi ' 

where the last step follows since < l in the last sum, with equality if and only 
if p k — for all k ^ i + 2. This implies m < oo and the uniqueness of p (given 
to) with the non-zero frequencies as claimed. □ 

It it possible to analyze the case of internal UC on the basis of the compact sets 
to be introduced below in Section|4] However, as J. Hofbauer pointed out to us 1 8 1, 
it is more natural to start with a larger compact set to be introduced in il Tj . Our 
main result in this section is thus 

Theorem 2. Assume that, for the initial condition p(0) and fixed r > 1, the r-th 
moment exists, X)/c>o k r Pk(®) < °°- Then, to = X)fc>o ^Pfc(O) infinite and, both 
in discrete and in continuous time, lim^oo \\p(t) — p\\ 1 — with the appropriate 
fixed point pfrom Proposition^ 

The proof relies on the following lemma, which slightly modifies and completes 
the convergence arguments of 1221 Sec. 4.1], puts them on rigorous grounds, and 
extends them to continuous time. 

Lemma 2. Let r > 1 be arbitrary, but fixed. Consider the set of probability mea- 
sures with fixed mean to < oo and a centered r-th moment bounded by C < oo, 

Mt, m ,c = (pe Mf : ]Tfcp fe = to, M r ( P ) s: C}, (11) 

equipped with (the metric induced by) the total variation norm, where 

M s (p) = Y.\ k - m \ s p* w 

k^O 

for s € This is a compact and convex space. Both M\ and M r satisfy 

M s (1Zq(p)) $J M s {p), with equality if and only if p is a fixed point of TZq- Fur- 
thermore, Mi is a continuous mapping from m c to K^o and a Lyapunov 
function for the dynamics in continuous time. 
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Proof. Let a sequence (p( n >) C M-^ m q be given and consider the random vari- 
ables = (fc) fcgNo on the probability spaces (No,p*- n - ) ). Their expectation 
values are equal to to, which, by Markov's inequality lETl p. 599], implies the 
tightness of the sequence (p^). Hence, by Prohorov's theorem |21 Thm. III. 2.1] 
(see also [3 Chs. 1.1 and 1.5]), it contains a convergent subsequence (recall 
that, by Proposition^ norm and pointwise convergence are equivalent in this case). 
Let p S A4i denote its limit and / = (k) k£fsSg a random variable on (No,p), to 
which the f( ni > converge in distribution. Since r > 1, the f( ni > are uniformly 
integrable by Markov's and Holder's inequalities. Hence, due to |9 Lemma 3.11], 
their expectation values, which all equal m, converge to the one of /, which is 
thus to as well. Consider now the random variables g( ni > = g = — m\ r ) keNo 
on (No,p'"') and (No,p), respectively. The expectation values of the g( ni *) are 
bounded by C, which, again by (9] Lemma 3.11], is then also an upper bound 
for the expectation value of g (to which the converge in distribution). This 
proves the compactness of Mt m c - The convexity is obvious. 
With respect to the second statement, consider 



\l — m\ 

PkPt 

I H- — M 

kAi^i^kVi 



kve 

PkPt o Z (\i~m\ s + \k + e~i~m\ s ). 



(13) 



- + \k-t\ 2 

k,£^0 1 1 i=k/\l 



For notational convenience, let j = k + I — i. We now show 

\i-m\ s + \k + £-i-m\ s < \k - m\ s + \l - m\ s . (14) 

If {k,£} = {i, j}, then HAi holds with equality. Otherwise, assume, without loss 
of generality, that k < i ^ j < I. If m ^ k or m ^ £, we have equality for ,s = 1 
but a strict inequality for ,s = r due to the convexity of x i— > x r . (For s = 1, this 
describes the fact that a recombination event between two sequences that are both 
longer or both shorter than the mean does not change their mean distance to the 
mean copy number.) In the remaining cases, the inequality is strict as well. Hence, 
M s (lZo(p)) ^ M s (p) with equality if and only if p is a fixed point of IZo, since 
otherwise the sum in (II 31 contains at least one term for which (114-i holds as a strict 
inequality. 

To see that Mi is continuous, select a converging sequence (p^) in A^^ m c 

and consider the random variables — (\k — m\) k£fs!g on (No,p^ ra - > ). As above, 
the latter are uniformly integrable, from which the continuity of Mi follows. Since 
M%(p) is linear in p and thus infinitely differentiable, so is the solution p(t) for 
every initial condition p e M-X m c> compare |1 Thm. 9.5 and Rem. 9.6(b)]. 
Therefore, we have 

Mi{p Q ) = liminf Ml(pW) 7 Ml(Po) - Mi(no(p )) - Mi(p ) < 0, 
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again with equality if and only if p is a fixed point. Thus, Mi is a Lyapunov 
function. □ 

Proof of Theorem^ By assumption, the r-th moment of p(0) exists, which is 
equivalent to the existence of the centered r-th moment by Minkowski's inequality 
1211 Sec. II.6.6]. This obviously implies the existence of the mean m. By Lemma|2l 
p(t) <E -Mi m c follows for all t ^ 0, directly for discrete time and via a satisfied 
subtangent condition 1 14, Thm. VI.2.1] (see also 1 1 Thm. 16.5]) for continuous 
time. In the discrete case, due to the compactness of A4f m c , there is a convergent 
subsequence (p(U)) with some limit p. Consider now the mean distance Mi to the 
mean copy number from H 2i . If lim^oo p(t) = p, we have, due to the continuity 
of Mi and TZq, 

MtCRoip)) = lim Mi(7Zo(p(*))) = lim Mi(p(t + 1)) = Mi(p), 

t — *oo t — *oo 

thus p is a fixed point by Lemma [2] Otherwise, there are two convergent sub- 
sequences (p(ti)), with limit p, and (p(si)), with limit q, whose indices alter- 
nate, ti < Si < ti + \. Then, we also have M\ (TZq (p(U))) ^ Mi(p(sj)) and 
Mi(n (p(si))) > Mi(p(t i+ i)), and therefore 

Afi(p) > M^Koip)) = lim Afi(72o(p(*i))) ^ lim Mi(p( Si )) = M x (g) 

Z— »00 i — ► OC 

> Mi(^ (q)) = lim Afi(Ho(p(«i))) ^ lim Mi(p(* i+1 )) = M x (j>). 

i — >oo i — >00 

Thus, both p and q are fixed points by Lemma[2]and hence equal by Proposition^ 
In continuous time, the claim follows from Theorem ^ since Mi is a Lyapunov 
function by Lemma|2] □ 

Note that, for q = 0, the recombinator can be expressed in terms of explicit 
frequencies irk t g of fragment pairs before concatenation (with copy numbers k and 
£) as TZo(p)i = Y^j=o ^jji-P However, we have, so far, not been able to use this 
for a simplification of the above treatment. 



4. Alternative probability representations 

Our next goal is to find the analogue of Theorem Efor the case of q = 1 (ran- 
dom UC). Whereas the convergence arguments for the case q — relied on a 
compact set of probability measure defined via the r-th moment, we are not (yet) 
able to extend this approach to q > 0. Instead, we will consider, as an alternative 
representation for a probability measure p € Mf, the generating function 

</Hz) = 5> fc z fc , (15) 

for which ^(1) = ||p||i = 1 and the radius of convergence is at least 1. We will 
restrict our discussion to such p for which lim sup^^ A/p^ < 1. Then, the radius 
of convergence, p(tp) = l/limsup^^ ^p^by Hadamard's formula 1 20] 10.5], 
is larger than 1. This is, biologically, no restriction since for any 'realistic' system 
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there are only finitely many non-zero p k (and therefore p(tp) = oo). Mathemati- 
cally, this condition ensures the existence of all moments and enables us to go back 
and forth between the probability measure and its generating function, even when 
looked at ip(z) only in the vicinity of z = 1 (see Proposition [6] below and 1211 
Sec. 11.12]). By abuse of notation, we define the induced recombinator for these 
generating functions as 

K{ij,){z) = Y,nP) k z k - (16) 

In general, with the exception of the case q — 1, we do not know any simple 
expression for lZ(ip) in terms of ip. Nevertheless, (I16> will be central to our further 
analysis. 

It is advantageous to use the local expansion around 2 = 1, written in the form 
^(z) = ^(fc + l)a fc (z-l) fc , (17) 

whose coefficients are given by 

* " (ITT)I' i, '" <1) " STlg Q * = : * »' <18 > 

In particular, a = 1 and a± = h J2e>o ^Pt This definition of a k is size biased, 
and will become clear from the simplified dynamics for q = 1 that results from it. 
For the sake of compact notation, we use a = (afc) fcgNo both for the coefficients 
and for the mapping. The coefficients a are elements of the following compact, 
convex metric space. 

Definition 3. For fixed a and S with < a ^ S < oo, let 

X a j = {a = (a k )keNo ■ «o = 1, oi = a, < a*, < S k fork > 2}. 
On this space, define the metric 

d(a,b) = Y, d k\ a k-h\ (19) 
with dk — (j/6) k for some < 7 < |. 

It is obvious that d is indeed a metric and that X a> s is a convex set, i.e., we have 
T) a + (1 — rf)b £ X at s for all a, b e X a ,s and 77 e [0, 1]. Note that we use the 
same symbol d as in Q since it will always be clear which metric is meant. The 
space X a ,s is naturally embedded in the Banach space (cf. 1261 Sec. 24.1]) 

H l/S = {xe R Nn : ||as|| < 00} (20) 

with the norm \\x\\ = 2^fc>o(T/^) fe l Xfe l' ^ or 7 an( l ^ as l n Definition|3] In partic- 
ular, d(a, b) = \\a — b\\. Furthermore, we have the following two propositions. 
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Proposition 5. The space X a ,s is compact in the metric d of dl9l . 

Proof. In metric spaces, compactness and sequential compactness are equivalent, 
compare 1 12, Thm. II.3.8]. Hence, let (aS n ') be any sequence in X a ,s. By assump- 
tion, Oq"' 1 = 1 and CLi = a. Furthermore, each element sequence (a^) C [0, S k ] 
has a convergent subsequence. We now inductively define, for every k, a conver- 
gent subsequence (a^ lk '^), with limit ctfc, such that the indices {n^.i : i £ N} 
are a subset of the preceding indices {nk-i,i '■ i £ N}. This way, we can pro- 
ceed to a 'diagonal' sequence (a < -" i •*>). The latter is now shown to converge to 
a = (flfc), which is obviously an element of X ai 5. To this end, let e > be 
given. Choose m large enough such that J2k>m(^) k < an d then * sucn that 

EfcLo 4l4 " a k\ < e/2. Then 

m 

d(a("-U) = E^i4 nM) - o*i + E d *i4 nM) 



(21) 



which proves the claim. □ 

Proposition 6. If limsup^^ Mp^ < 1, the coefficients from Jl 81 exist and 
a{p) £ X aj g with a = a(p) 1 = = | X)fc>o ^Pfe an< ^ some Conversely, if 
p(a) £ X a g for some a, S, one has limsup^^ typ^ < 1. 

For a proof, we need the following 

Lemma 3. Let fo(z) = X)fc>o ° kzk ^ e a P ower series with non-negative coeffi- 
cients Cfc and fx(z) = X)fc>o ~hJo^ ( x )( z ~ x ) k tne expansion °f fo around some 
x £ [0, p(fo)[. Then, p(fo) = X + p(f x ), including the case that both radii of 
convergence are infinite. 

Proof. Since the open disc B x (p{fo) — x) is entirely included in Bo(p(/o))> the 
inequality p(f x ) ^ p(fo) ~ x immediately follows from the theorem of rep- 
resentability by power series [2(5] Thm. 10.16]. Consider now the power series 

f xe * v (z) = Y, k >o vJo ] ( xelv )i z - xe%ip ) k with arbitrary <p £ [0,2w[. Its coef- 
ficients satisfy \ti k \xe^)\ < £ n ^ fe J^F)l c k xn ~ k = ffiix) due t0 the non- 
negativity of the Cfe. This implies p(f xe iv) ^ p(fx) by Hadamard's formula. 
Therefore, fo admits an analytic continuation on Bq(x + p(f x )), the uniqueness 
of which follows from the monodromy theorem [20 Thm. 16.16]. The theorem of 
representability by power series then implies the inequality p(fo) x + p(f x ), 
which, together with the opposite inequality above, proves the claim. □ 

Proof of Proposition]^ The assumption implies p(i/j) > 1 for ip from J15I . Then, 
from Lemma[3] we know that limsupj,^^ y/{k + l)afc < oo. Since furthermore 
a-k ^ (k + l)afc, also limsup^^ tyak < oo, so there is an upper bound 6 
for ^/ok and thus a(p) £ X aj g. The converse statement follows from Jl 8i and 
Lemma[3] □ 
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Therefore, any mapping from X a j into itself that is continuous with respect to the 
metric d from dl9l has a fixed point by the Leray-Schauder-Tychonov theorem 
EH Thm. V.19]. 

Note further that each X a ^ contains a maximal element with respect to the 
partial order a ^ b defined by a k ^ b k for all k £ No, which is given by 
(1, a, 8 2 , (S 3 , . . .). This property finally leads to 

Proposition 7. The space P a ^ := {p £ : a(p) £ X a ,s}, equipped with (the 
metric induced by) the total variation norm, is compact and convex. 

The proof is based on the following two lemmas. 

Lemma 4. For any subset of P a ,8, the corresponding generating functions from 
(I15i are locally bounded on i? 1+1 /,j(0). 

Proof. It is sufficient to show boundedness on every compact K C Bi + i/s(Q), see 
llT9l Sec. 7.1]. Thus, let such a Kbe given and fix r £ [0, 4[ so that K is contained 
in £>i +r (0). Then, for every p £ P a s and every z £ K, 

= |X>* z *| ^ E^ 1 + r ) fc = ^ + r ) 

= ^(fc + l)a(p) fc r fe < l + 2ar + ^(fc + l)(n5) fe < 00, 

where r<5 < 1 was used. This needed to be shown. □ 

Lemma 5. If for a sequence (p^) C P a ,s, the coefficients = a(p( n ' > ) from 
dl8> converge to some a with respect to the metric dfrom d!9i . then the generating 
functions ip n from H5\ converge compactly to some ijj with ip(z) = Y2k>oPk zk 
and the p^ thus converge in norm to p £ P a ,s- 

Proof. By to Lemma|4] the sequence (ip n ) is locally bounded in S 1+1 /i (0). Due 
to the pointwise convergence \a k n ^ — a k \ ^ d k 1 d(a^ n \ a) — > 0, we have 

^ fe )(l) = (fc + l)!ai" ) ^ (fc + l)!a fe = ^)(l) 

for every k £ No- Then, the compact convergence ip n — ► i/j follows from Vitali's 
theorem 1191 Thm. 7.3.2]. In particular, this implies that — ■> p k ^ and 

Now, choose r G ] 1, 1 + | [■ Then there is, for every e > 0, an n e such that 
sup| z |^ r \^{z) — 4> n {z)\ < £ for all n ^ n e . This implies 



\P k n) -P k \ 



— I M^m dz 



for all n ^ n s by Cauchy's integral formula 1 13, Thm. 7.3]. Now, let e > be 
given. Then 

\\p {n) -p\\ 1 = Y,\p^-P^< £ T L i 

fcJsO r 

for all n > n e , which proves the claim. □ 
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Proof of Proposition [7] Let (p^) denote an arbitrary sequence in P a j and 
(o") = (ofpW)) the corresponding sequence in X a t $. Due to Proposition [5] 
there is a convergent subsequence (a' n '") i . Then, by Lemma|5] (p' 71 *^ converges 
in norm to some p E P a .s- This proves the compactness property. The convexity 
of P a} s is a simple consequence of the convexity of A4f, the linearity of the map- 
ping a, and the convexity of X a □ 

Another property of the mapping a : P a j — > X at $ is stated in 

Lemma 6. For every a and 8, the mapping a : P a j — > X a g from (II 81 is contin- 
uous (with respect to the total variation norm and the metric d) and injective. Its 
inverse p : a(P a j) — > P a j is continuous as well. 

Proof. Let p,q <G P a .s and assume a(p) = a(q). Then, as in the proof of 
Lemma |3 the uniqueness of the generating function in f? 1+1 /,5(0) follows, and 
thus p = q, which proves the injectivity of a. The other statements follow from 
Vitali's theorem 1191 Thm. 7.3.2]: Norm convergence of a sequence (p^ n ') in 
P a: s to some p implies convergence of its element sequences and thus compact 
convergence of the corresponding generating functions ip n to ip, which is given 
by ip(z) — Ylk>oPk z ■ This, in turn, implies convergence of each sequence 
(a(p^)k) to a(p) k , from which, as in i21\ . the convergence (a(p(™))) — > a(p) 
(with respect to d) follows. The converse is the statement of Lemma |5] (see also 
EU Prop. 1.6.8]). □ 

Note that, if p(ip) > 2, the inverse of the mapping a is given by 
p(a) k =J2(-iy- k ( e k )(t+l)a e . 



5. Random unequal crossover 

Let us now turn to the random UC model, described by q = 1 in 0. Here, the 
recombinator simplifies to 1221 (3.1)] 

l + min{k,£,i,k + £-i} 

n ^ = E o — {k + m+1) — w 

k+l^i 

As for internal UC, by Lemma[0 the reversibility condition JlOb directly leads to 
an expression for fixed points, 



Pk Pi _ Pi Pj 



k+l£+l i+lj+1 



for all k + £ = i + j . 



This has p k = C(k + l)x k as a solution, with appropriate parameter x and nor- 
malization constant C. Again, it turns out that all fixed points are given this way. 
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Proposition 8. |22 Thm. A. 2] Every fixed point p g of 7Z\ is of the form 



where m = J2k^s0 kp k ^ 0. □ 

One can verify this in several ways, one being a direct inductive calculation. 
The main result of this section is 



Theorem 3. Assume that lim sup^.^^ \/p fc (0) < 1. Then, both in discrete and in 
continuous time, limt^oo \\p(t) — p\\i — with the appropriate fixed point pfrom 
Proposition^ 

For a proof, we consider the following alternative process, verbally described in 
11221 p. 720f]. It is a two-step stick breaking and glueing procedure which ulti- 
mately induces the same (deterministic) dynamics as random UC, even though the 
underlying process is rather different. This will lead to a simple expression for the 
induced recombinator of the coefficients a from d 1 St t, which allows for an explicit 
solution. 

Proposition 9. Let p € M.f. Then, 

^ = J2ih Pe (24) 

gives a probability measure 7r G and the recombinator can be written as 

i 

Ki (p)i = 77 j = ( n * n )i > ( 25 ) 

3=0 

where * denotes the convolution in ^(Nq). 

Here, ( 1241 describes a breaking process in which, without any pairing, each se- 
quence is cut equally likely between any two of its building blocks. In a sec- 
ond step, described by d25i . these fragments are paired randomly and joined (or 
'glued'). 

Proof. It is easily verified that n is normalized to 1. With respect to ( 1251 . note the 
following identity for k + I ^ i, 

\{j : [i -I) V < j ^ i A k}\ = 1 + min{fc, £, i, k + i - i}, 

which can be shown by treating the four cases on the LHS separately. With this, 
inserting (|24j into the RHS of (j25j yields 

' 1 

X>i E {k + m+1) p*p< 

3=0 j=0 k^j t^i-j y n ' 
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^ iAfe 

= E ( k + i)(e + i) PkPe E 1 

El + min{fc,£, i,k + i — i} 

This nice structure has an analogue on the level of the generating functions. 

Proposition 10. Under the assumptions of Theorem^ let 4>{z) = X]fc>o ^k 2 ^ 
denote the generating function for it from i24l . Then, 

1 f 1 

^ Z ) = — Z J and R iW(«) = ^) - 

Proof. Recall that ip(z) = J2k>oPk zk - Equations (1241 and (125 > lead to 

= E E ZTI^ = £ ZTT^ E - E 7Ti^ : T^T~ 

k^oe^k e^o k^e e^o 

and, due to absolute convergence of the series involved, 

ft x (v>)(*) = n ^p)kZ k = E z " E 

= E ^ E n k-ez k ~* = <i>{z? ■ 

i^O k^l □ 

The following lemma states that the radius of convergence of ip does not decrease 
under the random UC dynamics. Thus, it is ensured that, if p(ip) > 1, also lZi(ip) 
may be described by an expansion at z = 1, i.e., by coefficients a. 

Lemma 7. The radius of convergence of TZi(ip) is p(lZi(?p)) p(tp)- 



Proof. As l/p(ip) = limsup^^ >fp^ =: x < 1 and limfc_oo y/c+T = 1, 
there is a constant C > with p fe ^ C(fc + for all k. Note the identity 

n 

+ min{i, J, n — i,n — j}) = (i + l)(n — i + 1) 



for i ^ n, which follows from an elementary calculation. Then, 1221 implies 

^i(p)*^C 2 ^ x fe+£ (l + min{fc,£,i,/c + ^-i}) 
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C 2 E x " E*- 1 + min{i,j,n — i,n — j}) 
3=0 

2 



1 1 ^ i 



C 2 (i+l)f^+l)x'= (- 



C 



(i + l)x i 



Accordingly, limsupj,^^ y/Hi(p) k ^ x ^ 1, which proves the claim. □ 

These results enable us to derive the following expression for the coefficients 
a, using the expansion of illi : 



z-1 

k>0 



<H0 dC 

2 



(26) 



= E E a « 

fe>0 \n=0 



Ofc- 



So it is natural to define the induced recombinator 



1 fc 



(27) 



n=0 



for which we have 



Lemma 8. 77ze recombinator TZ\ given by Mil maps each space X a j into itself 
and is continuous with respect to the metric dfrom ( 119k 

Proof. Let a, <5 > be given and a,b £ X a> g, Trivially, TZi(a) = 1 and 
72.1(a)! = a. For k ^ 2, TZi{a) k = ■ k -^iJ2e^k a e a k-e ^ $ k ■ Tnis proves the 
first statement. For the continuity, note first that every TZ\{a) k with k ^ 2 is con- 
tinuous as a mapping from X a< s to [0, 5 k ] . Now, let e > be given. Choose n large 
enough so that J2k>n C^l) k < e /2> where 7 is the parameter introduced in Defini- 
tion!!] Then, there is an 77 > such that ^=2(1 1 5 ) k \Ki{a)k ~ ^i( b )fc| < e/ 2 
for a,b £ X a j with d(a, b) < 17. Thus, for such a and b, 

d(^(o),^i(6)) < e (J)" i^( a ) fe -^iWfci + E( 2 ^) fc < £ > 



fc=0 



k>n 



which proves the claim. 



□ 



Note that the fixed point equation on the level of the coefficients a is always satis- 
fied for ao and a±. If k > 1, one obtains the recursion 



^ ft— 1 

Ofe = T J" E anak -" 



which shows that at most one fixed point with given mean can exist. 
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Let us now consider the case of discrete time first. Analogously to (jSJi, define 
a(t) = a(p(t)) as the coefficients belonging to p(t), which are assumed to exist. 
It is clear from (O, (|26} and (27) that a(t + 1) = Ki(a(t)). We then have the 
following two propositions. 

Proposition 11. Assume a(0) to exist. Then, in discrete time, 
lim aj.it) — a for all k 0. 

t — >oo 

This result indicates that a weaker condition than the one of Theorem may be 
sufficient for convergence of pit). 

Proof. Clearly, ao(t) = 1, ai{t) = a. Furthermore, by the assumption and (1261 . 
the coefficients <Zfc(t) exist for all k, t S No- Now, assume that the claim holds for 
all k ^ n with some n and let k = n + 1. According to the general properties of 
lim sup and lim inf , we then have 

1 k 

limsup a,k(t+ 1) < \imswp(ae(t)ak-e (t)) 



k 

-a 



.. — limsup Ofc(t) 

1 ft + 1 t^oo 

and analogously with for lim inf. This leads to 

^ limsup a k (t) < < |^limmfa fe (i), 

from which the claim follows for all k ^ n + 1 and, by induction over n, for all 
jfe 0. □ 

Proposition 12. 77;e recombinator TZ\, acting on X a $, is a strict contraction with 
respect to the metric dfrom \\9\ , i.e., there is a C < 1 such that, for all elements 
a,b e X a j, 

d(Hi(a),ni(b)) ^ Cd(a,b). 

Proof. First consider, for k ^ 2, without using the special choice of the d k , 

k 



k^2 e=o 
1 I k 

= j_ | 1 |y^( a ^ _ bi)(a k -£ + bk-i) 



k^2 i=0 

k>2 K + 1 1=2 

l>2 k>l 1 



(28) 
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With the choice d k = (j/5) h , where we had 7 < i, we can find, for I ^ 2, an 
upper bound for the inner sum, 

k>e k>i 1 



which, together with J28L proves the claim. □ 

Together with Banach's fixed point theorem (compare 1181 Thm. V.18]), the two 
propositions imply that a(t) converges to (1, a, a 2 , . . .) with respect to the metric 
d, and that convergence is exponentially fast. 

In continuous time, we consider the time derivative of a(t) :— a(p(f )), which 
is, by l[T8}, 

f t a(t) = f t a(p(t)) = a(fti (p(t)) - p(t)) = 7^ (a(t)) - a(t) . (29) 

The following lemma ensures, together with [ 1 Thm. 7.6 and Rem. 7.10(b)], that 
this initial value problem has a unique solution for all a(0) = ao 6 X a ,5. 

Lemma 9. Consider the Banach space H^/g from J20L vv/f/i some < 7 < 4, 

and its open subset Y = {x G -ff-y/a : |a3fc| < (25) }. 77ien, the recombinator IZi 
from \21\ maps Y into itself satisfies a global Lipschitz condition, and is bounded 
on Y. Furthermore, it is infinitely differentiable, IZi € C°°(Y, Y). 

Proof. For x £ Y, one has |xfc| < (2<5) fc , hence |72.i(a;) fc | < (25) k , with a similar 
argument as in the proof of Lemma|H] Consequently, 72. 1 (Y) C Y. So, let x, y E Y. 
Then, similarly to the proof of Proposition^] one shows the Lipschitz condition 

\\Hi(x) - Ki(y)\\ < Q'l*' ~ Wl E fcTT^^ < T^2" I|:E ~ y|1 
and, since < 1/(1 — 27) in Y, the boundedness, 

ll*x(«)||< T ^IW|<^^ 5 . 
With respect to differentiability, consider, for sufficiently small h € Y, 

2 

7ti(a; + h) k = Ui{x) k + — ^i t - ( h e + Hi(h) k . 



Since 



it is clear that 7£.i is differentiable with linear (and thus continuous) derivative, 
whose Jacobi matrix is explicitly it' 1 (x) kl = (%rHi(x) k — j^xu-t if k I 

and zero otherwise, hence one has TZi € C 1 (F, F). It is now trivial to show that 
Hi € C 2 (y, y) with constant second derivative and thus Hi <G C°° (Y, Y). □ 
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Proposition 13. //ao £ X a ^$ for some a, 5, then a(t) € X a ^$ for all t ^ and 

lim^oo d(a(t), a) = with a — (1, a, a 2 . a 3 , . . .). 

Proof. The first statement follows from [ 14 Thm. VI.2.1] (see also [l] Thm. 16.5]) 
since, due to the convexity of X Q < 5, we have a + t(lZ(a) — a) € X Q $ for every 
a € 5 and t € [0, 1], hence a subtangent condition is satisfied. For the second, 
observe that IZi (ct) = ex. We now show that 

L(a ) = d(a ,a) (30) 

is a Lyapunov function, cf. Definition^ With the notation of Lemma[9] note that 
the compact metric space X ai s is contained in the open subset Y of the Banach 
space H~jg. The continuity of L is obvious. Now, let ao € X a< s be given. By 
Lemma|5|and [ 1 Thm. 9.5 and Rem. 9.6(b)], the solution a(t) of d29i is infinitely 
differentiable. Thus, for t E [0, 1], 

L(a(t)) - L(a Q ) = ||a + t(Ki(a ) - a ) + o(t) - a\\ - ||a - a\\ 

< t(||^i(a )-^i(a)|| - ||a +o(t), 

where o(t) is the usual Landau symbol and represents some function that vanishes 
faster than t as t — > 0. From this, by the strict contraction property of IZi (Propo- 
sition ! 12> . the Lyapunov property l|8} follows, with equality if and only if ao = a. 
Since X a t g is compact, TheoremHimplies the claim. □ 

We are now able to give the previously postponed 

Proof of Theorem^ By Proposition [6j we have a(0) = a(p(0)) £ X a> s with 
a = |m and some 5. In discrete time, according to Propositions ^2 ar, d H3 
and Banach's fixed point theorem (compare 1 18 Thm. V.18]), it then follows that 
a(t) — * a = (1, a, a 2 , . . .) with respect to the metric d Inserting d23l > into Jl 8i 
and letting x — m/(m + 2) yields 

*-g ( ,-»w»W '-»> , (*+'>''' 

The claim now follows from Lemma [3] Similarly, in continuous time, the claim 
follows from Proposition^] □ 

Let us finally note 

Proposition 14. For the dynamics described by \29\ , the fixed point a from Propo- 
sition W^ is exponentially stable. 

Proof. Let ao € X a j be arbitrary. The Lyapunov function from the proof of 
Proposition^]satisfies, as a consequence of (13 1 1 and Proposition! 121 

L(a ) < d(lZi(a ), TZi(a)) — d(a , a) ^ — (1 — C) d(a , a), 

with < C < 1. From this, together with (|30j and (TJ Thm. 18.7], the claim 
follows. □ 
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Remark. In a related UC model introduced by Takahata 1 24 1, for which 

1 



T, 



ij,k£ — Oi+j,k+l 



k + e + i' 



the recombinator IZi appears for the coefficients b(p) k = (k + 1) a(p) k , where 
b(p) 1 is the mean copy number m. The above results then imply, under the appro- 
priate condition on p(0), that b(t) — > (1, to, to 2 , . . .) as t — > oo both in discrete 
and in continuous time. This corresponds to convergence of p(t) to the fixed point 

P with Pfc = ^TT(^TT) fc - 

6. The intermediate parameter regime 

In this section, q may take any value in [0,1]. With respect to reversibility of fixed 
points, one finds 

Proposition 15. For parameter values q £ ]0, 1[, any fixed point p G A4f of the 
recombinator lZ q , given by and Q, satisfies p k > for all k ^ (unless it 
is the trivial fixed point p = (1, 0,0,.. .) we excluded). None of these extra fixed 
points is reversible. 

Proof. Let a non-trivial fixed point p be given and choose any n > with p n > 0. 
Observe that T^ ln _ 1 m > for < q < 1 and hence 

Pn±l = ^j(p)n±l = T^i^ktPkPi ^ ^n+ln-l.nnPnPn > 0- 

The first statement follows now by induction. 

For the second statement, evaluate the reversibility condition dlOt for all com- 
binations of i, j, fc, I with i+j = fc + i?^4. This leads to four independent 
equations. Three of them can be transformed to the recursion 

(fc + l)<? p x ,,„,,, 
Pk= 2(k-l)+2 q p- Pk - 1 > ke ^^> 

from which one derives explicit equations for all p k with k e {2, 3, 4} in terms 
of p Q and p x . Inserting the one for p 2 into the remaining equation yields another 
equation for p 4 in terms of p and p lt which contradicts the first equation for all 
q G ]0, 1[, as is easily verified. □ 

So, non-trivial fixed points for < q < 1 are not reversible, and thus much more 
difficult to determine. Our most general result so far is 

Theorem 4. If p(0) G P a ,sfor some a, S, then p(t) G P a .s for all times t G No, 
respectively t € K^o> Ttq na $ a fixed point in P a ,8- 



The proof is based on the fact that lZ q is, in a certain sense, monotonic in the 
parameter q. This is stated in 
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Proposition 16. Assume a(p) S X a ,s for some a, S. Then, with respect to the 
partial order introduced before Proposition^ a(lZ q (p)) ^ a(lZ q >(p)) for all 
^ q ^ q' ^ 1. In particular, a(lZ q (p)) S X a .$ for all ^ q ^ 1. 

To show this, we need three rather technical lemmas. The first one collects formal 
conditions on the difference of two discrete probability distributions T^ ke with 
different parameter values (but j = k + I — i and the same fixed k, £), These are 
then verified in our case. 

Lemma 10. Let the numbers X{ € R (0 ^ i $J r with some r € No) satisfy the 
following three conditions: 

r 

^Xi = 0. (32) 

i=0 

x r -i — Xi for all ^ i ^ r. (33) 

{Xi ^ : O^i^n 
, i ■ (34) 
Xi < : n < i < [§J 

Further, let /i 6 1 (0 ^ i ^ r) fee given vw'f/z 

0</i-/o</a-/i<...</r-/r-i. 05) 

TTzen, we /zave 

r 
i=0 

Proof. Let us first consider the trivial cases. If Xi = 0, everything is clear, so 
let Xi ^ 0. If r ^ 1 then Xi = 0, so let r ^ 2, and thus n ^ § — 1. Define 
ii = ft = for odd r. Then, we can write 

i—0 i—0 i—n-\-l 

Furthermore, for r — i i, due to ( I35> . 

/i + /r-j = + /r-i+1 + (/i — /j-l) — (/r-i+1 — /r-i) ^ + /r-i+1 ■ 

Now, define C := 2~_^™=o ^ = — $^[=1+1 x i ~ > 0> an( ^ th e c l a i m follows 
with (1361 . since r — n^n + lby assumption: 



fi%i ^ C [/n ~t~ fr — n fn+1 fi — n— l] 



i=0 



= C [(/ P _„ - / r -„-l) - (/„+! - /„)] > 0. □ 



Lemma 11. Let j € No be fixed and fi — (£),., i £ No, where is the falling 
factorial, which equals 1 for j — and i(i — 1) • • ■ (i — j + 1) for j > 0, hence 
nzijT] f or i 3 s j- Then condition (I35t z'i satisfied. 
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Proof. For j = 0, condition (I35> is trivially true. Otherwise, each /; is a poly- 
nomial of degree j in i with zeros {0, 1, . . . , j — 1}, hence we have the equality 

= ft — fo = ... = fj-t — fj-2- Then, for i ^ j — 1, the polynomial and 
all its derivatives are increasing functions since lim^oo f = oo. Therefore, for 

1 > j — 1, we have ^ /j+i — ^ .^+2 — /i+i- Hence d35i holds. □ 

Lemma 12. For ^ q ^ q' ^ 1 and all k,£, Equations J32l i-l l34t are true for 
r = k + £and x, = T$ - T$, where T$ = T$ H with j = k + £-i. 

Proof. The validity of i32\ and J33i is clear from the normalization and the 
symmetry of the T^} . For J34> . let k ^ £ without loss of generality. In the trivial 

cases q = q' or k = 0, choose n — [§J- Otherwise, xi — T^J — T^j < for 

k ^ i ^ L5J ' smce Cm < , and xo > 0. For ^ i ^ fc, consider 

u- = — + 1 = I - 

1 ikl vy 

Here, the first factor is less than 1, the second is equal to 1 for k = i, greater than 
1 for ^ k < i, and strictly decreasing with i. Since x\ ^ if and only if yi ^ 1, 
there is an index n with the properties needed. □ 

Proof of Proposition^^ We assume ^ q ^ q' ^ 1, Lemmas ITTTM 1 21 imply, for 

all k,£,j G N with k + £ ^ j, 



/ , (i-jy.ikl iki 

Then, since T^j = for i > fe + £, 

= O'+l)! X/ E W7)\ T ikl 

< (7TT)T E ^ E (l4r^ } = W)i ' 

From this, together with Lemma|H] the claim follows. □ 

Proof of Theorem^ According to Proposition [HO lZ q maps P a> g into itself, and 
thus, in discrete time, p(t) € P a ,d for every i G No- The analogous statement is 
true for continuous time t € K^o- To see this, consider P aj $ as a closed subset of 
I 1 . Recall that lZ q — 1 is globally Lipschitz on i 1 by Proposition^ Moreover, for 
any p 6 P a s and t € [0, 1], Proposition0tells us that 

p + t(ft,(p) - p) = (1 - t)p + «J,(p) € P Q , 5 ■ 

This implies the positive invariance of P a j by 1141 Thm. VI.2.1] (see also Q] 
Thm. 16.5]). The existence of a fixed point once again follows from the Leray- 
Schauder-Tychonov theorem 1 18 Thm. V.19]. □ 
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On the basis of the above analysis, and further numerical work done to inves- 
tigate the fixed point properties 1171221 . it is plausible that, given the mean copy 
number m, never more than one fixed point for 7Z q exists. Due to the global con- 
vergence results at q = and q = 1, any non-uniqueness in the vicinity of these 
parameter values could only come from a bifurcation, not from an independent 
source. Numerical investigations indicate that no bifurcation is present, but this 
needs to be analyzed further. 

Furthermore, the Lipschitz constant for 7Z q can be expected to be continuous in 
the parameter q, hence to remain strictly less than 1 on the sets X a j in a neighbor- 
hood of q = 1. So, at least locally, the contraction property should be preserved. 
Nevertheless, we do not expand on this here since it seems possible to use a rather 
different approach |7|, which has been used for similar problems in game theory, 
to establish a slightly weaker type of convergence result for all < q < 1, and 
probably even on the larger compact set A4^ m c of Lemma[2] 

7. Concluding remarks 

In this article, we have shown that, for the extreme parameter values q = (in- 
ternal UC) and q = 1 (random UC), any initial configuration satisfying a specific 
condition converges to one of the known fixed points, both in discrete and contin- 
uous time. The condition to be met is, for q = 0, the existence of the r-th moment 
(r > 1, see Theorem|2j, respectively, for q = 1, that the corresponding generating 
function has a radius of convergence p > 1 (Theorem[3j. Convergence takes place 
in the total variation norm in all cases. As argued in the previous section, similar 
results can be expected for the intermediate parameter values as well. 

These results are valid for deterministic dynamics and thus correspond to the 
case of infinite populations. With respect to biological relevance, however, we add 
some arguments that it is reasonable to expect this to be a good description for 
large but finite populations as well, i.e., for the underlying (multitype) branching 
process. For finite state spaces, such as in the mutation-selection models discussed 
in |6|, the results by Ethier and Kurtz |4 Thm. 11.2.1] and the generalization |2 
Thm. V.7.2] of the Kesten-Stigum theorem 1 10 11 1 guarantee that in the infinite 
population limit the relative genotype frequencies of the branching process con- 
verge almost surely to the deterministic solution (if the population does not go 
to extinction). Since for the UC models considered here the equilibrium distribu- 
tions are exponentially small for large copy numbers (owing to Theorem |4] also 
for q £ ]0, 1[), one can expect these systems to behave very much like ones with 
finitely many genotypes. This is also supported by several simulations. Neverthe- 
less, this questions deserves further attention. 
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